Proton Spin Structure in the Resonance Region 
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We have examined the spin structure of the proton in the region of the nucleon resonances 
(1.085 GeV < W < 1.910 GeV) at an average four momentum transfer of Q 2 = 1.3 GeV 2 . Using the 
Jefferson Lab polarized electron beam, a spectrometer, and a polarized solid target, we measured 
the asymmetries A\\ and A± to high precision, and extracted the asymmetries A± and A2, and the 
spin structure functions g± and 32- We found a notably non-zero A±, significant contributions from 
higher-twist effects, and only weak support for polarized quark-hadron duality. 



PACS numbers: 13.40.Gp, 13.88.+e, 14.20.Dh 

Ever since the first polarized EMC experiment found 
that the proton's spin is not fully carried by its valence 
quarks [1], the nucleon spin structure has been studied 
extensively, for example at SLAC [2], CERN [3], and 
DESY [4]. The main focus has been on kinematics in 
the deep inelastic scattering (DIS) region and with large 
momentum transfer (high Q 2 ) where observations can be 
readily interpreted in a perturbative QCD framework. In 
more recent years, lower energy regimes (Q 2 ~ 1 GeV 2 ) 
have grown in importance, where the transition from the 
asymptotically free to the bound configuration of the 
quarks can be probed. Evaluating the requisite moments 
of the spin structure functions requires data or model de- 
pendent extrapolations, up to parton momentum fraction 
1=1, where the electromagnetic scattering probe inter- 
acts with the proton as a whole, including the kinematic 
region dominated by nucleon resonances. 



The new experimental focus on the larger x and lower 
Q 2 regimes so far has been concentrated on longitudi- 
nal polarization as the dominant, and technically more 
accessible, component. This has limited investigations 
to the spin structure function g\, which represents the 
charge-weighted quark helicity distributions. Due to the 
success of DIS interpretations based on the naive par- 
ton model, which is limited to longitudinal spin com- 
ponents, little emphasis has been placed on transverse 
spin studies. In contrast, the operator product expan- 
sion (OPE) approach to QCD includes transverse spin 
starting from leading twist [5, 6] and is applicable at 
all kinematics. Since each higher order of twist, inter- 
pretable as increased correlation between partons, adds 
another l/Q term, higher twist contributions should be 
more prominent at low Q 2 . Although transverse spin is 
suppressed in DIS at leading twist, thus motivating the 
twist- 2 WW approximation [7] 
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9 r(x,Q 2 ) = -gi(x,Q 2 )+ I 9l (y,Q 2 )^, (1) 



twist-3 should contribute significantly to g 2 . The OPE 
relates the twist-3 matrix element d 2l which represents 
quark-gluon correlations, to the third moments of g\ and 
92 [6], 

d 2 = 3 f x 2 (g 2 -g 2 vw )dx = [ x 2 {2 9l + 3g 2 ) dx, (2) 
Jo Jo 

permitting comparison with QCD lattice calculations, 
and thus providing a clean test of the theory. 

Also, the phenomenon of quark-hadron duality [8] has 
captured much interest [9, 10, 11, 12, 13]. This concept 
connects the pcrturbative QCD description of quarks 
and gluons to the hadronic description at lower energies. 
Whereas global duality relates the entire resonance re- 
gion to extrapolations of the asymptotic structure func- 
tions, local duality relates the extrapolated quantities to 
restricted resonance regions. Both have been observed in 
unpolarized scattering [14], but precision measurements 
of the spin structure functions in the resonance region are 
needed to determine whether they also display duality. If 
local duality, as first seen by Bloom and Gilman [8] , were 
observed for g\ , this would not only demonstrate its uni- 
versal, rather than accidental nature [15], but it could 
also serve to justify extrapolations. The potential useful- 
ness of the latter for the determination of the moments 
increases at low Q 2 , where the resonances extend over 
wider ranges of x than at high Q 2 . 

Experiment E01-006 was conducted in Hall C of the 
Thomas Jefferson National Accelerator Facility by the 
Resonance Spin Structure (RSS) collaboration. Utilizing 
established procedures and equipment, we have measured 
the asymmetries An and A± in the scattering of polarized 
electrons off a polarized proton target. These asymme- 
tries are defined as the dimensionless, relative difference 
between the cross sections with parallel and anti-parallel 
(or perpendicular and anti-perpendicular) alignment of 
the proton and the electron spins. Focusing exclusively 
on the hadronic vertex, we are left with the virtual pho- 
ton asymmetries A\ and A 2 , which are functions of the 
virtual photon's four-momentum-squared — Q 2 and the 
invariant mass W of the final state. Without a measure- 
ment of A± , this step would require the use of a model or 
a sweeping assumption. Further, switching from an ex- 
ternal scattering view to an internal structure interpreta- 
tion, we can obtain the spin structure functions gi{x, Q 2 ) 
and g 2 (x,Q 2 ), where Q 2 is interpretable as the energy 
scale set by the probe and, in DIS, x represents the frac- 
tion of the proton momentum carried by the constituent 
the probe interacted with. The exact relation between 
these quantities, and other relevant definitions, can be 
found in Sec. II of Ref. [2] or in Sec. 2.1 of Ref. [9]. 



The experiment used Jefferson Lab's continuous, po- 
larized electron beam with energy of 5.755 GeV and a 
nominal current of 100 nA. The beam polarization was 
measured by a M0ller polarimeter [16] installed upstream 
of the target, and the beam hclicity was flipped at 30 Hz 
on a pseudo-random basis. 

Frozen 15 NH3, in 1 — 2 mm fragments, was used as 
the proton target in the University of Virginia appara- 
tus [17] in which a 4 He evaporation refrigerator at 1 K 
coupled with a 5 T polarizing magnet created a stable 
polarization environment. The polarization population 
enhancement was achieved via microwave pumping (dy- 
namic nuclear polarization) and measured by an NMR 
system using pickup coils embedded in the target mate- 
rial. For the A± measurement, the entire target appa- 
ratus was rotated by 90° in the scattering plane. About 
equal amounts of data were taken with each polarization 
direction, flipping the nuclear spins by adjusting the mi- 
crowave frequency, to reduce systematic effects. To main- 
tain uniform polarization in the bulk target material, the 
beam was continually moved across the face of the target 
in a 1 cm maximum radius spiral raster pattern around 
the beam axis. This extra degree of freedom required a 
dedicated beam position monitor [18] for accurate event 
reconstruction. 

Scattered electrons were detected using the High Mo- 
mentum Spectrometer (HMS), positioned at a scattering 
angle of 13.15°. Two different HMS momentum settings 
were used, 4.078 and 4.703 GeV, to cover the desired wide 
range m W, resulting in 0.8 < Q 2 < 1.4 GeV 2 . A detec- 
tor package consisting of hodoscope planes, wire cham- 
bers, a gas Cerenkov counter, and a lead glass calorime- 
ter allowed for particle identification and measurement 
of the event kinematics. A more detailed description of 
the apparatus and technique can be found in Ref. [19]. 

Approximately 160 million scattering events were 
recorded on the proton target, resulting in highly precise 
determinations of the parallel and perpendicular asym- 
metries. These are obtained from observed raw event 
counting asymmetries which are scaled to 100% polar- 
ization and corrected for contamination from radiative 
and dilution processes: 
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Here, N^ 1 is the charge corrected observed yield 
for the parallel (perpendicular) and anti-parallel (anti- 
perpendicular) spin alignment, respectively. Pf, and Pt 
are the beam and target polarizations, / is the dilution 
factor, Cn is a small 15 N nuclear polarization correction, 
and fnc and Arc are radiative corrections. 

The corrected asymmetries A\\ and Aj_ are shown in 
Fig. 1; Aj_ is notably non-zero. The average proton po- 
larization was 62% (70%), and the beam polarization was 
71% (66%) during the parallel (perpendicular) running. 
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FIG. 1: Our measured asymmetries A\\ and A±, fully cor- 
rected (points) and without radiative corrections (curves). 



For the parallel alignment, the product Pi, x P t was de- 
rived by normalizing the measured elastic asymmetry [19] 
to the known value, resulting in better accuracy than 
achievable from direct measurements. In the perpendic- 
ular case, the limited knowledge of the elastic asymme- 
try made the direct measurement of and Pt the better 
choice. The systematic errors are summarized in Table I, 
highlighting the lack of models and data for perpendicu- 
lar radiative corrections. 

The dilution factor represents the fraction of events 
that truly scattered from a polarized proton in the tar- 
get. It was determined from the ratio of free proton to 
total target rates calculated via a Monte Carlo simula- 
tion which had been matched to calibration data acquired 
specifically for this purpose. The QFS parameterization 
[20], modified to improve agreement with our data, was 
used as input for the Born inelastic cross sections for 
A > 3 nuclei. Fits to Hall C inelastic e-p data [21] were 
used for the H contribution. The unpolarizcd structure 
function F\ and the ratio R of the longitudinal to trans- 
verse cross sections are derived from the same fits [21, 22]. 
The uncertainty in these cross sections was the dominant 
source of systematic error for the dilution correction. 

Convoluting radiative prescriptions with models of the 
resonance region, the elastic peak, and our target, we 
obtained radiated cross sections and asymmetries. The 
external radiative corrections were determined using the 
procedure established in [23], while the POLRAD soft- 
ware [24] was used to determine the internal radiative 



TABLE I: Averaged systematic errors in the asymmetries. 



Error Source 


Ay 


A ± 


Target Polarization 
Beam Polarization 
Dilution Factor 
Radiative Corrections 
Kinematic Reconstruction 


}1.1 % 

4.9 % 
2.7 % 
0.4 % 


2.9 % 
1.3 % 
4.9 % 
12.9 % 
0.4 % 



corrections. The resonance fit model was iteratively im- 
proved, until the radiated values matched our experimen- 
tal data. The model then trivially provided the correc- 
tions to our measurements, with fnc accounting for the 
radiative dilution from the elastic tail and Arc for all 
other influences. 

We extract the virtual photon asymmetries A\ and A2, 
shown in Fig. 2, from the corrected physics asymmetries 
Am and A±, using only R as model dependent input. 
The spin structure functions g\ and gi (Figs. 3 and 4) 
are then obtained utilizing F\. The uncertainties in F\ 
and R are included in our total systematic error and in 
the error bands of our plots. 

We have fitted the W dependence of our A\ and A2 
data using an approach similar to that applied to un- 
polarized cross sections in Ref. [23], substituting for the 
DIS component a form based on the phenomenological 
spin structure paramctcrizations of Ref s. [25, 26]. These 
fits served as input in the iterative procedure to obtain 
our radiative corrections and to calculate the integrals of 
<?i and <?2 at constant Q 2 . Each spin asymmetry was fit- 
ted independently, since they represent different physical 
quantities. 

To test quantitatively for global duality in g\, we can 
integrate in x over the resonance region and compare 
the results obtained from resonance data and DIS ex- 
trapolations (Fig. 3). We used our fit to integrate over 
the resonance region (1.09 < W < 1.91 GeV) and took 
the average of the integrals from several DIS extrapola- 
tions calculated from target mass corrected [27] , next-to- 
lcading order parton distribution functions (NLO PDFs) 
[28, 29, 30] over the same range of W at our aver- 
age Q 2 = 1.3 GeV 2 . We found the ratio of integrals, 
PDFs to data, to be 1.17 ±0.08, indicating agreement at 
only the two sigma level, and suggesting that PDF ex- 
trapolations into the resonance region may not be valid 
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FIG. 2: Virtual photon asymmetries A\ and A2 from our 
data and corresponding fits. Also shown is the E155 fit to 
DIS data [25, 26], evaluated at our (x,Q 2 ), and the Soffer 
limit for A2 [31], based on our Ai fit. The upper error band 
indicates the systematic error in Ai, the lower one Ai- 
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FIG. 4: Our (RSS) values for gi and the approximation g^ 
(Eq. 1) as evaluated from our data. 



FIG. 3: Results for g± from this experiment (RSS) and other 
relevant data [2, 32], as well as target mass corrected NLO 
PDFs. The upper scale shows W (at Q 2 = 1.3 GeV 2 ) for 
reference. 



at this Q 2 . The ratios for restricted, but still rather 
broad, W ranges differ from unity by several sigmas, 
demonstrating that local polarized duality is not valid 
at our Q 2 : 1.09 < W < 1.4 GeV is 6.47 ± 0.95, and 
1.4 < W < 1.91 GeV is 0.87 ± 0.06. Including also 
large x resummations for the PDFs [11], the global ra- 
tio changed to 1.42 ± 0.10. The quoted errors are based 
on the data integrals only, including a 0.4% contribution 
from computing our fit at fixed Q 2 . Our results are in 
good agreement with the recent results from CLAS [13]. 

Approximate global duality within errors was reported 
by [10], based on A\ resonance data averaged over a 
broad Q 2 range from 1.6 to 2.9 GeV 2 and compared to a 
DIS fit to data. The weak Q 2 dependence of A\ (within 
large errors) allows for averaging, instead of calculating 
the ratio at each Q 2 value as is required for testing du- 
ality in the structure functions. But duality in the spin 
asymmetry A\ oc gi/F\ could be due to accidental can- 
cellations in the ratio gi/F\. 

Our results for g2 are much clearer, especially in the 
framework of the QCD OPE. The comparison of our data 
and the g™ w approximation, evaluated from our mea- 
surements of gi , provides strong evidence of the signifi- 
cance of higher-twist terms at this Q 2 , as shown in Fig. 4. 
Combining our measurements of g\ and g^ we can inves- 
tigate specifically the twist-3 contribution via the matrix 
element d 2 (Eq. 2). Over the measured range (0.29 <x< 
0.84), we find ~d 2 = 0.0057 ±0.0009 (stat) ±0.0007 (syst), 
including a 4% contribution to the systematic error from 
our fit's assumed Q 2 dependence. This significantly non- 
zero result highlights the limitation of leading-twist ap- 
proximations. Extrapolating this result to Q 2 = 5 GeV 2 , 
assuming a 1/Q dependence, we find d 2 = 0.0029 com- 
pared to the SLAC result d 2 = 0.0032 ± 0.0017 [33]. 

In summary, our results significantly increase the avail- 
able information on the proton spin structure: These new 



data provide a connection to the measurements at DIS 
kinematics and fill a significant void in the explored re- 
gions. Our measurement with transverse spin arrange- 
ment is the first in the resonance region, with notably 
non-zero results. Our data clearly indicate the impor- 
tance of higher twist contributions and thus quark-gluon 
correlations. We have established that Bloom-Gilman 
polarized duality is meaningful only for the resonance re- 
gion as a whole, although local polarized duality may yet 
be observed at higher Q 2 ranges. 
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APPENDIX 

The results of our fits to A\ and A2 are detailed in 
table II and plotted against our data in figure 2. They 
are based on the method of Ref. [23] , combining a num- 
ber of Breit-Wigner shaped resonances BWi and a deep- 
inelastic tail: 
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We have optimized the values of amplitude a*, centroid 
Wj, and width gi of four Breit-Wigner resonances, as well 
as the exponent a and the polynomial coefficients f3 n for 
the DIS tail. The values of the parameters Xi , k , and ji 
were kept as in Ref. [23] . While the fits are a function of 
W, the DIS contribution is evaluated at the mean values 
of x and Q 2 of our data at that W bin. 

In addition to the explicit contribution in the terms Ki 
and K cm , and the indicated 1 / y/Q 2 term for the DIS part 
of A2 , our fit function depends on Q 2 via the interrelation 
of W, x and Q 2 : 



TABLE II: Results of Fits to Ai and A 2 , each based on 28 
data points as plotted in figure 2. 



Parameter 


Ai Fit 


A 2 Fit 


x7dof 


1.24 


1.42 


ai 


-0.562 ± 0.133 


-0.147 ± 0.052 


a 2 


0.401 ± 0.069 


0.220 ± 0.049 


03 


0.537 ± 0.042 


0.158 ± 0.029 


a4 


0.246 ± 0.044 


0.173 ± 0.029 


toi 


1.204 ± 0.010 


1.214 ± 0.021 


IV2 


1.345 ± 0.005 


1.338 ± 0.010 


W3 


1.544 ± 0.012 


1.479 ± 0.021 


W4, 


1.734 ± 0.009 


1.653 ± 0.015 


01 


0.161 ± 0.065 


0.137 ± 0.086 


S2 


0.063 ± 0.025 


0.071 ± 0.035 


93 


0.283 ± 0.058 


0.113 ± 0.082 


34 


0.125 ± 0.043 


0.129 ± 0.062 


Of 


0.031 ± 0.391 


0.458 ± 0.375 


jflb 


0.186 ± 0.074 


0.100 ± 0.040 


/3i 


-0.032 ± 0.158 


0.094 ± 0.079 


/3 2 


-0.393 ± 0.272 


-0.119 ± 0.133 




0.957 ± 0.403 


-0.025 ± 0.202 



The choice of four resonances resulted in better fits 
than using only three, an option also supported by ex- 
amining the second derivatives of the data smoothed with 
a spline curve. Additionally, the fits as reported here also 
gave excellent agreement with the cross sections cttt' and 
<7lt' obtained from our data and a ctt model [22] . 

Table III details the results of duality integral ratios 
for various sub-ranges (those matching the resonances 
assumes in our fit to A\, as well as those of [9]) and 
for the different PDFs [28, 29, 30]. The PDF integrals 
over individual resonances differed significantly from the 
data value, indicating that at our Q 2 , local duality is not 
observed in g\. 



TABLE III: Ratio of NLO PDFs to data, integrated over the 
indicated ranges. The top five W ranges are based on the 
location of the Breit-Wigner shapes resulting from our fit to 
Ai and the last two match those of Ref. [9] ; all are bounded 
by our acceptance. The elastic contribution [9] was calculated 
using the dipole form of the form factors. 





W Range 


BSB 


GRSV 


AAC 


Average 


Delta 


1.09-1.30GeV 


3.41 


4.18 


3.96 


3.93T0.58 


Rl 


1.30-1.39GeV 


1.28 


1.44 


1.33 


1.36T0.10 


R2 


1.39-1.68GeV 


0.77 


0.82 


0.75 


0.78±0.05 


R3 


1.68-1. 79GeV 


0.77 


0.84 


0.77 


0.79T0.06 


Global 


1.09-1.91GeV 


1.11 


1.23 


1.14 


1.17T0.08 


Delta + Rl 


1.09-1.40GeV 


5.78 


6.84 


6.44 


6.47T0.95 


R2 + 


1.40-1.91GeV 


0.84 


0.91 


0.84 


0.87±0.06 



W 2 = M 2 - Q 2 + 



The resonance terms are otherwise (^-independent. 



